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^ ■ ABSTRACT 

Aims. We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS 
survey and the physical properties of their host galaxies. We study a bolometric regime (<Lboi > 8 x 10^^ erg s~^) where several 
theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. 

Methods. To derive robust estimates of the host galaxy properties, we use an SED fitting technique to distinguish the AGN and host 
galaxy emission. We evaluate the eff'ect on galaxy properties estimates of being unable to remove the nuclear emission from the SED. 
The superb multi-wavelength coverage of the COSMOS field allows us to obtain reliable estimates of the total stellar masses and 
star formation rates of the hosts. We supplement this information with a morphological analysis of the ACS/HST images, optical 
spectroscopy, and an X-ray spectral analysis. 

Results. We confirm that obscured quasars mainly reside in massive galaxies (M^ > 10^^ M©) and that the fraction of galaxies hosting 
such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude 
diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. 
We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical 
properties of the hosts. We find that at z~ 1, ^ 62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable 
to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ^ 71% at z~ 2, and 
100% at z~ 3. We also find that the evolution from z~ 1 to z~ 3 of the specific star-formation rate of the Type-2 QSO hosts is in 
excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the 
objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated 
galaxies tend to host Type-2 QSOs with low Eddington ratios (A < 0.1), while disk-dominated or merging galaxies have at their 
centers BHs accreting at high Eddington ratios (A > 0.1). 



Key words. (Galaxies:) quasars: general - Galaxies: nuclei - Galaxies: active - Galaxies: evolution - Galaxies: star formation - 
X-rays: general 



1. Introduction 

There is strong observational evidence that the formation 
and growth of supermassive black holes (SMBHs) and 
their host galaxies are related processes: e.g. the tight 
corr elation between the SMBH mass and the luniinos- 
itv (iKormendv & Richsto"iie] 1 19951: 'McL ure & Dunlop I l20Q2l : 
iMarconi&HuntI l2Q03h. mass (lM agorriai 7et al. I Il998h. ve- 
locity dispersion (Ferrarese & Merritt 2000; "Gebhardt et alTI 
12000 : Merritt & Ferrarese 2001; Tremaine et al. 2002), or the 
structure (iGraham et al. 1 12001) of the bulgeQ. Additional evi- 
dence of the connection between the two phenomena is pro- 
vided by the similarity between the global star-formation rate 
(SFR) and active galactic nuclei (AGN) activity: both peak 
at 2 and decline rapidly at lower ( e.g. Dickinson et al 



20031: iMerloni et al. 1 12004) and higher (e.g. lWilkins et al. 11200 



Brusa et al. redshifts. From the theoretical side, a num 

ber of models have been proposed to link the formation and 
evolution of SMBHs to the structure formation over cosmic 
time (e.g .Kauffmann & Haehn elt 2000; Somerville et al. 2001 



Granato et al. 1 120041: iMonaco & Fontanot 2005; Granato et al. 



20061: LCroton et al."^ '20061: iLapietal.i 12006: .Hopkins et al 
20061: [Cen 2011). Some of thes e semi-analytical mod els and 



hydrodynamical simulations (e.g. ' Springel et al.1l2005b invoke 
major mergers of gas-rich galaxies as the mechanism en- 
abling the fueling of the central SMBHs and the building 
of the galaxy's bulge. Alternative fueling mechanisms have 
been discussed in the literature, including minor-mergers (e.g. 
[Johansson et al. I l2009l). bars ( e.g. IJogee I I2004D . disk instabil- 
ities (e.g. iGenzel et al. I l2008l) . and recycled gas from dying 



^ See'Peng ' (120071) andUa hnke & Maccio I (1201 Ih for alternative ideas 
on the origin of these scaling relations. 
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Fig. 1. X-ray luminosity in the [0.5-10 keV] band corrected for 
obscuration plotted against the Nh column density for the X-ray 
sources in the XMM-COSMOS survey with a redshift. The 142 
Type-2 QSOs studied in this paper are indicated by filled circles. 
The dashed lines delimit the area where L[o.5-iO)tey] > 10"^"^ erg 
s"^ andNn > 10^^ cm'l 



stars (e.g. ICiotti et al. Il2010h . iHopkins & Hemquist I (l2009h ar- 
gued that such alternative processes may be enough to fuel 
Seyfert-like AGNs, but that extreme mechanisms such as ma- 
jor mergers are necessary to provide the gas supply needed for 
bright quasars. Further more, some observational studies in- 
dicate that there is a morphological transition in AGN hosts 
around Lboi ~ lO^^L© between a mixture of disk and ellipti- 
cal galaxies to exclusively bulge-dominated galaxies (see Fig. 
1 of iHopkins & Hemquist I t2009 and references therein). In this 
framework, the interest in studying the host galaxy properties of 
bright quasars is clear. These o bjects also dominate the space 
density of AGN at 2 (e.g. U eda et al. Il2003l: iHasinger etaTI 
2005). An obvious complication in the study of the host galaxy 
properties of quasars is usually the emission of the central 
AGN which outshines the galaxy light making it extremely dif- 
ficult to derive constraints on the colors, stellar populations, 
and morphologies of the host. The situation clearly improves 
if the powerful quasar is obscured by gas and dust on parsec 
scales: we refer to these objects as Type-2 QSOs. Although the 
fraction of obscured AGN is found to decrease with lumi nos- 
ity in several studies (X-ray: [Ueda et al. 2003, Steff'en e t al. I 
L2003", 'Hasin2er""200< Hasinger 2008, Brusa et al. 2010; op- 
tical: Simpson 2005: MIR: Maio lino et al. 2007 : Treister et^alTI 
l2008h . a non-negligible population of obscured QSO s is still re- 
quired by X-ray background synthesis models (e.g. 'Gilli e t alTI 
2007). Radio-loud Type-2 QSOs h ave been o bserved for many 
years thanks to radio surveys (see iMcCarthy Ill993l for a com- 
prehensive review), while radio-quiet Type-2 QSOs were de- 
tecte d more recently in Chandra and ^KM^-lSlewton X-ray sur- 
veys (Dawson et al. '2001; "Norman et al. 2002; Mainieri e t al. 
2002; Stern et al. 2002; Delia Ceca et al. 2003; Fiore et al. 



200 31: IPerola et al. II2004I: iTozzi et al. 2006; Barger et al. 2005; 
iMateos et al. I l2005l: iKmmpe et al? ^2008; VignaH et al ...201Q) 
optica l surveys (SDSS, |Zakamska et al. 2003, Re ves et al. I 
20081). and mid -IR observations ( Martinez- Sansigr e et al. H2QQ6[ 



Polletta et aO2 008: Lanzuisi et al. 2009). The aim of this paper 



is to extend to quasar luminosities the study of the host galaxy 
properties of AGN, using a large sample of X-ray selected Type- 
2 QSOs. In Section 2, we describe the selection criteria, while 
in Sections 3 and 4 we present the average spectral properties of 
the sample in the X-ray and optical band, respectively. Section 5 
describes the method used to separate the host galaxy from the 
AGN emission. In Section 6, we illustrate the derived properties 
for the host galaxies, and in Section 7 discuss their morphologi- 
cal appearance. Finally, in Section 8 we summarize the results of 
this work. We quote in this paper magnitu des in the AB sy stem, 
we use a Chabrier Initial Mass Function (IChabrier 1120031) . and 
we assume a cosmology with Hq = 70 km s"^ Mpc"\ Qm = 0.3, 
and Qa = 0.7. 



2. Sample selection 

We started from the sample of ~ 1 800 p oint-like X-ray sources 
detected in the XMM-COSMOS survey (ICappelluti et al. II2009I) 
for which optical identifications, multi wavelength photometry, 
and a compilation of available redshifts from ongoing spec- 
troscopic campaigns in the COSMOS field (Lilly et al. 2 009; 
Trump et al. 2009) were presented in iBrusa et aP ('2010'). For 
the remaining sou rces very a ccurate photometric redshifts are 
available ( Sal vato et al. ]12009'). For each X -ray source, we per- 
formed a detailed X-ray spectral analysis to quantify their X-ray 
luminosity and obscuration. The procedure adopted to extract 
source and background spectra is equivalent to that described 
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Fig. 2. Redshift distribution of the Type-2 QSO sample. The 
hatched histogram is for spectroscopic redshifts, while the empty 
one includes both spectroscopic and photometric redshifts. 



by iMainieri et al. 1 (l2007h . We used XSPEC0 (vl 1.3.2) to per- 
form our spectral fitting analysis. We first fitted the EPIC-pn 
data with two basic input models: a single powerlaw (PL) and 
a powerlaw modified by intrinsic absorption at the redshift of 
the source (APL). The PL model consists of two XSPEC compo- 
nents wabs'^zpowerlw, while the APL model consists of the com- 
bination of three diff'erent components wabs^zwabs^zpowerlw. 
In both models, the wabs component accounts for the photo- 
electric absorption due to the Galactic column density and is 
fixed to the average value in the COSMOS region (N^^^ ~ 
2.7 X 10^^ cm"^ iDickev & Lockman1ll99Qlfl . The zwabs com- 
ponent describes th e photoelectric absorp t ion us ing Wisconsin 
cross-sections (Morrison & McCammon L Il983b and the only 
two parameters are the equivalent hydrogen column density and 
the source redshift. zpowerlw is a single power-law parametrized 
by the photon index, the redshift, and a normalization factorQ 
Additional components are included to describe the presence of 
soft-excess or the Fe Ka emission line. The model fits yield the 
power-law photon index F, and the X-ray luminosity in both the 
[0.5-2] and [2-10] keV rest- frame bands, and the APL model also 
provides the intrinsic column density Nh. 

We used the following selection criteria to select Type-2 



QSOs: a) de-absorbed Lx [0.5 - 10 keV 
Nh > 10^^ cm~^ (see also IMainieri et al. 



> ly^ erg s-^ b) 
2002t ISzokolv et al. I 



l2004l:lTozzi et al. ll2006l: lKrumpe et al. "2008"). We excluded four 
sources (XID= 82, 215, 261 8, 5208) that would formally match 
these criteria but for which iBrusa et aiTI (1201 Ol) proved that the 
XMM-A^^w^6>^ emission is the blend of two Chandra sources. 
The final sample of X-ray selected Type-2 QSOs consists of 142 
object^. 

In Fig.[Tl these candidates are plotted in the X-ray luminos- 
ity versus column density plane. The redshift distribution of the 
Type-2 QSO sample is shown in Fig. O 33 (^ 23%) of them 
have a spectroscopic redshift. For this sub-sample of sources, we 
verified that no broad lines (FWHM> 2000 km s"^) are present 



^ http : //heas arc .gsfc.nasa. go v/doc s/xanadu/xspec/ 
^ This is an average value for the Galactic Nh in the COSMOS area, 




channel energy (keV) 




where N^""^ is in the range [2.5-2.9] xlO^^ cm~^. This range in Galactic 
column density does not affect the results of our spectral analysis. 
^ We refer the reader to http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ ^^1^^ 



channel energy (keV) 

Fig. 3. Stacked X-ray spectrum of the 33 Type-2 QSOs with 
a spectroscopic redshift (left) and 79 Type-2 QSOs with only 
a photometric redshift (right). Adjacent bins were combined 
until they had a significant detection at least as large as 
1.5(T. The solid line is the best fit to the continuum with a 
wabs^zwabs^powerlaw XSPEC model. In the left panel, clear 
residuals are present around the rest-frame 6.4 keV energy. 



in the optical spectrum. For the remaining 77% of the sample, 
we u sed the very accurate photometric redshifts ( ISalvato et al. I 
[200^ . The X-ray selected sample of Type-2 QSOs presented in 
this work is complementary to the SDSS opticall y selec t ed one 
presented by Zakams ka et al. . (.2003 ) and Reyes et al. I (|2008|) 
in terms of redshift range covered: the SDSS sample based on 
[OIII] luminosities cannot be extended to redshifts higher than 
0.8, while our sample is almost entirely at z> 0.8. 
From the absorption-corrected luminosities in the rest-frame 
[2-10 keV] energy band, we estimated the bolometric luminosi- 
ties for the Type-2 QSOs using the l uminosity-dependent bolo- 
metric correction factor (fe-ioy^ev) of Marconi et al. I (|2004|) . For 
the X-ray luminosities sampled by our Type-2 QSOs, k2-iokeV 
is in the range [25,120], and the bolometric luminosities cover 
the range Uoi = [1.2 x ^"^^4.5 x 10"^^] erg s"^ with a mean 



<Lboi >~ 8 X 10^^ erg s . A limitation of this 



for further details on the spectral models 

5 



Source XlD-2028 discussed in Brusa e t aTl j201Qh is not included 
in this paper because a re-analysis of the X-ray spectrum returned a 
column density value of log (Nh)= 21.7 + 0.1 cm~^. 



approach is the intrinsic spread in AGN SEDs, which could 
lead to significant diff'erences in the bolometric correction from 
that predicted based on a mean energy distributi on for AGN 
(see lVasudevan & Fabian 1 120071: iLusso et al. iflOld) . From these 
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bolometric luminosities, and assuming a standard radiative effi- 
ciency of 6rad = 0.1, we derive a median accretion rate onto the 
black holes, M = Lboi/^rad^^^, of approximately ^ 1 Mo/year, 
simi lar to what observed for optically selected Type-1 QSOs 
(e.g. iKoUmeier et ani2Q06l: iHopkins et ani2Q09l) . 

3. Stacked X-ray spectrum 

In order to study the average X-ray properties of our Type-2 
QSO sample, we obtain a stacked X-ray spectrum using the 
EPIC-pn data. To avoid the possibility of the blurring and broad- 
ening of the Fe Ka line caused by the uncertainties in the pho- 
tometric redshifts, we first considered only the 33 sources with 
a spectroscopic redshift. The background- subtracted data of in- 
dividual objects in the energy range corresponding to the rest- 
frame 2-10 keV were divided into 40 intervals so that each in- 
terval corresponds to a rest-frame 200 eV. These spectra were 
corrected for the detector response and summed together. The 
details of the procedure will be presented in a future paper 
(Iwasawa et al. in preparation). The smallest net count recorded 
in the rest- frame 2-10 keV for the spectroscopic sample is 28 
photons. The resulting stacked spectrum (2682 net counts) is 
shown in the left panel of Fig. [3] It was further rebinned for 
display purposes: adjacent bins were combined until they had a 
significant detection at least as large as l.Scr. From the stacked 
spectrum, we measured a column density Nh = ^-^I'l ^ 1^^^ 
cm"^ (quoted errors are Icr), which is consistent with the aver- 
age value of <Nh >= 5.4 x 10^^ cm"^ obtained from the fits to 
the single X-ray spectra. As is clear from the left panel of Fig. 
[21 an absorbed power-law model leaves large residuals around 
the expected energy of the Fe Ka line. If we add a gaussian 
model centered on 6.4 keV in the XSPEC fitting procedure, the 
quality of the fit improves significantly with 99% confidence ac- 
cording to an F-test. For the Fe Ka equivalent width, we mea- 
sured a value of EW= 104 eV. Taking into account the known 
anti-correlation between the Fe Ka and the X-ray luminosity 
(llwasawa & Taniguchilll993h . this value i s consistent with typ- 
ical EW values for Co mpton-thin AGNs (iBianchi et al. 1 120071: 
IChaudharv et"ani2010 '). 

We repeated this procedure for the photometric redshift sam- 
ple selecting objects that have more than 30 counts in the rest- 
frame 2-10 keV band. This selection leaves 79 objects out of 
112 and the corresponding stacked X-ray spectrum (6,005 net 
counts) is shown in the right panel of Fig. [3] Owing to the 
larger number of objects, the final spectrum has more counts and 
consequently the uncertainties in the column density values are 
smaller, Nr = ^-^I'l x 10^^ cm"^. At the same time, the signifi- 
cance of Fe Ka is significantly lower, and we ascribed this to the 
uncertainties, although small, in the photometric redshifts that 
tend to blur out the signal. In fact, an uncertainty in the photo- 
metric redshift value as low as = 0.02x(l-i-z) would move the 
centroid of the 6.4 keV line outside the rest- frame energy bin of 
200 eV used in Fig.S 

4. Optical spectroscopy 

Of the 33 Type-2 QSOs with a spectroscopic redshift, 31 have 
an optical spectrum from the zCOSMOS survey, while the re- 
maining two were observed with the IMACS spectrograph at 
Magellan. We used the 31 zCOSMOS spectra to construct an av- 
erage optical spectrum for the Type-2 QSO sample. We grouped 
the spectra into redshift bins to ensure that each bin had a com- 
mon spectral range. The VIMOS low-resolution (LR) spectra ob- 




1000 2000 3000 4000 

Wavelength (angstroms) 



Fig. 4. Stacked optical spectrum of 31 Type-2 QSOs from the 
zCOSMOS survey. The LR and MR spectra have been grouped 
in redshift bins, as specified in the legend inside the figure, in 
order to have a common wavelength coverage between them. 
The more evident emission/absorption lines are labeled in the 
figure. 



tained with the LR-blue grism (3600-6800 A) were divided into 
two redshift bins: seven objects in 2.4 < z < 3.2 and four objects 
in 0.8 <z< 1.5. An additional twenty medium-resolution (MR) 
spectra (5550-9650 A) covering the redshift range 0.6 < z < 1.6 
were stacked together. The stacked optical spectrum covers the 
rest-frame wavelength range from 900 to 4500 A (see Fig. |4). 
Both [NeV]3426 and [Oil] 3727 are clearly detected, and the ra- 
tio of the two lines is [OII]/[NeV]~ 4.1. A similar value for 
the lines ratio was measu red for a [OIII150Q7 sel ected sample 
of SPSS Type-2 Q SOs bv lZakamska etal.l (l2003l) . In addition, 
iGim et al. I (l2010l) found a ratio [OII]/[NeV]> 4 for most of the 
obscured quasars selected from the SDSS using the ratio of the 
2-10 keV flux to the [NeV]3426 emission line flux. 
We considered two features in the stacked optical spectrum of 
Fig.m the 4000 A break, D„(4000), and the absorption line. 
Because of the wavelengths of these two features, the follow- 
ing discussion is limited to the 20 Type-2 QSOs in the redshift 
range 0.6 < z < 1.6 with a VIMOS-MR spectrum, hence the 
mean optical spectra may not be representative of the overall 
population of Type-2 QSOs. The D„(4000) is a rough stellar age 
indicator, with small values of the break for young stellar popu- 
lations and large values of the index for old ste llar populations. 
Following the definition of lBalogh etan(ll999h . we derived the 
p„(4000) as the ratio of the fluxes measured in the [3850-3950] 
A to [4000-4100] A bands. In the first band, we subtracted from 
the stacked spectrum the emission associated with the [Nelll] 
line. The second age indicator is based on the absorption line, 
i.e. the H^a index defined bv lWorthev & Ottaviani! (Il997h : the 
continuum is estimated from the average flux in the two bands 
[4041.60,4079.75] A and [4128.50,4161.00] A, and the index 
is expressed in terms of an equivalent width integrating within 
the band [4083.50,4122.25] A. A strong absorption line is 
a signature of a large fraction of young stars, hence a recent 
burst of star formation. We measured D„(4000)= 1.19 ± 0.02 
an d H^A = 4.7 ± 0.4. We can compare these values with Fig. 7 
of IWinter et al. I (l2010l) . These authors used the two spectral fea- 
tures as age indicators in the study of local AGN selected in the 
14-195 keV band with SWIFT. Our values fall in a region where 
the spectral models have a significant (> 30%) contribution from 
a young stellar populations. 
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Fig. 5. Two examples of SED decomposition. Black circles are 
the observed photometry in the rest- frame. Purple and blue lines 
correspond respectively to the galaxy and AGN template found 
as the best fit solution, while the black line shows their sum. 



5. Distinguishing the AGN and host galaxy 
emission with SED fitting 

The main goal of our study is to use the rich multiwavelength 
coverage of the COSMOS field to constrain the properties of the 
galaxies hosting Type-2 QSOs. In particular, we wish to esti- 
mate their stellar mass and star-formation rate (SFR). The tech- 
nique we used is a detailed model fitting of the total SED of 
the Type-2 QSOs for 14 diff'erent photometric bands that en- 
compass optical to MIR wavelengths: six SUBARU bands (B, 
V, g, r, i, z); U, J, and K bands from CFHT; four Spitzer/IRAC 
channels, and lAjdm from Spitzer/MIPS (see iBrusa et al. |[2QTqI 



and references therein). This allows us to sample a wide wave- 
length interval ranging from ~ 3 800 A (Ucfht) to 24 jim. As al- 
ready presented in Merloni et al. (2010), we fitted the observed 
SED with a grid of models compiled by combining AGN and 
host gal axy templa t es. Fo r the AGN component, we adopted 
the Rich ards et aTl (l2006h mean QSO SED derived from the 
study of 259 IR selected quasars with both Sloan Digital Sky 
Survey and Spitzer photometry. We allow for extinction of the 
nuclear AGN light by applying a SMC-like dust-reddening law 
(iPrevot etal. 1 11984') of the form: AaIE{B -V)= 1.39 /l;^^. The 
E{B - V) values for the AGN component are allowed to vary in 
the range 1 < E{B - V) < 9. Assuming a Galactic gas-to-dust 
ratio, our selection criteria Nr > 10^^ cm"^ w ould correspond 
to E{B - V) > 1.7, but there is evidence (e.g. iMaiolino et al. I 
2001) that the E(B-V)/Nh ratio in the circumnuclear region of 
an AGN may be lower than the Galactic value, at least for part 
of the AGN population. Our choice to set the lower boundary 
to E{B - y) = 1 is a compromise between allowing lower than 
Galactic values for the dust-to-gas ratio and avoid degeneracies 
in the SED fitting procedure caused by the inclusion of an unob- 
scured AGN template. Spectropolarimetry observations of Type- 
2 QSOs from the SDSS (Zakamska et al. , 2005) and of narrow 
line radio galaxies ( Tadhunter e t al. « .2002) show that in some 
sources there could be a non-negligible constribution from scat- 
tered quasar light. Assuming th at the s cattering process is in- 
dependent of wavelength ( Vernet et al. I, 2001), the SED of an 
unobscured quasars should be a good representation of this addi- 
tional component. To quantify the eff'ect of AGN scattered light 
on the estimates of the host galaxies properties presented in the 
next section, we repeated the SED fitting with < E{B - V) < 9 
for the AGN component. The values for the host stellar mass are 
essentially unchanged: for 97% of the objects the diff'erence in 
stellar mass is less tha n 5%. Stellar popula tion synthesis (SPS) 
models, as the one of i Bruzual & Chariot I (^003) used in this 
paper, combine individual simple stellar population spectra to 
match the observed broad band photometry of a galaxy. Old stel- 
lar populations are the main contributor to the total galaxy stel- 
lar mass, and their SED dominates at wavelength above ^ 600 
nm. This explains why adding the SED of an unobscured AGN, 
that would contribute significantly mostly in the UV, does not 
aff'ect the stellar mass estimates based on SPS models. As for 
the star-formation rate, we find that scattered AGN light would 
change significantly the values of SFR for 15% of the objects. 
Therefore, in spite of the possible multi-component nature of the 
optical continuum, the SED fits unambiguously require the pres- 
ence of a young stellar population for the majority, ^ 85%, of the 
objects. We also verified that the spread on SFR values, gener- 
ated from the inclusion or not of an unobscured AGN component 
in the SED fitting procudure, does not aff'ect the overall proper- 
ties of the host galaxies population (see Fig. 8,9, 10, 1 1). As for 
the upper boundary in the reddening for the AGN component, 
E{B - y) < 9, it corresponds to the average X-ray column den- 
sity of our Type-2 QSO sample of Nh ~ 5 x 10^^ cm"^ for a 
Galactic dust-to-gas ratio. We verified that allowing higher val- 
ues of E(B-V) would not aff'ect the values for the physical quan- 
tities derived from the SED fitting, since an AGN template with 
E{B- y) = 9 has pr actically no significant impact below ^ 2yum. 
While lMerloni et al. I (|2010), dealing with unobscured QSOs, as- 
sumed that all the 24 iim flux is due to the AGN, we tried to es- 
timate t he AGN contr i jbution at those wavelengths using the re- 
sults of iGandhi et al. I (12009^ . These authors found a strong cor- 
relation between the core luminosity at 12.3 jim and the X-ray 
luminosity in the [2-10 keV] band for a sample of local Seyfert 
galaxies. In the mid-infrared, they used high-resolution obser- 
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vations with the VISIR/VLT instrument, and owing to the su- 
perb angular resolution of those images were able to minimize 
the contamination by unresolved extended emission to the core 
fluxes. We used their Eq. (2) for the best-fit correlation obtained 
by considering only the 22 well-resolved sources 

iogi^^Y^) = (0.19 ± 0.05) + (1.11 ± om)log(^^^^). 

Using the formula above, we estimated the AGN flux at 123 
fim from the L[2-10 keV] luminosities corrected for obscuration 
obtained from the X-ray spectral analysis. We associated an er- 
ror in this mi d- infra red flux using the dispersion observed by 
iGandhi et al. I (l2QQ9h in the ratio of mid-infrared to X-ray lumi- 
nosity, cr(logLMiR/Lx) = 0.18 (see their Table 2 and Fig. 2). 
Finally, we included in the SED fitting procedure the constraint 
that the AGN component matches this predicted flux in the mid- 
infrared. 

For the host galaxy, we created a library of synthetic spec- 
tra using_Jhe_well^^ population synthesis models 
of iBruzual & Chariot] (l2003h . We compiled ten exponentially 
declining star-formation histories (SFH) SFR oc ^"^age/r ^j^j^ 
e-folding times (r) ranging from 0.1 to 30 Gyr, plus a model 
with constant star formation. For each of these SFH, we calcu- 
lated the synthetic spectrum at difl'erent ages, ^age, ranging from 
50 Myr to 5 Gyr, subject only to the constraint that the age 
should be younger than the age of the Universe at the redshift 
of the source. Finally, we included dust exti nction, modeled us- 
ing the Calzetti's law (ICalzetti et al. I[2000h. with values in the 
range < E(B - V) < 0.5. FoUowingTonta na et al. I (12004 ) and 
iPozzetti et aD (l2007h . we imposed the prior E(B - V) < 0.15 if 
kg^/^ > 4, a significant extinction being allow ed only for galax- 
ies with a high SFR. We adopted a Chabrierl (|2003„) initial mass 
function (IMFfl 

Finally, an additional component that could contribut e to the UV 
conti nuum of the SED is nebular continuum ( see Dick son et al. I 
Il995l) . For a sample of te n narrow line radio gal axies, possibly 
radio-loud Type-2 QSQs. iTadhunter et al. I (l2002l) calculated the 
fractional contribution of this componenent to the UV flux in 
the range 3 - 30%. In this work, we will not consider the possi- 
ble contribution of the nebular continuum. We note that at least 
for stellar mass estimates the efl'ect of this component should be 
negligible: the intensity of the nebular continuum drops quickly 
with wavelength and becomes extremelly faint at A > 500nm, 
where we expect to be more sensitive to the galaxy stellar mass. 

For each one of the Type-2 QSOs, we derived the best fit 
SED from the minimization of the comparing the observed 
and template fluxes at the redshift of the source. The combina- 
tion of templates allowed us to disentangle the emission of the 
host galaxy from the contribution of the central black hole. We 
show in Fig.[5]two examples of this composite fit. From the best 
fit SED, we used the galaxy component to estimate both the stel- 
lar mass and SFR of the host. For each one of these physical 
quantities, we estimated the confidence region corresponding to 
68% of confidence level for one in teresting parameter, impos- 
ing thatx^ ~ ^Ibest-fit - ^ (lAvni 111976. and Table 1 therein). 
The best-fit values and errors obtained are reported in Table 1 . 
We flagged in the same table the sources without a reliable fit to 
their SED for various reasons (e.g. photometry contaminated by 
a nearby bright object). 

^ Different choices of IMF lead to systematic shifts in the estimated 
stella r masses for any given SED, with the maximum shift for a Salpeter 
IMF dSalpeter 1955) given by M^^saipeter ^ M^^chabrier x 1-7 (see e.g. 
Pozzetti et al. 2007; Ilbert et al. 2010). 
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Fig. 6. Comparison between the estimates of (top panel) and 
SFR (bottom panel) obtained with the SED fitting technique us- 
ing only galaxies templates (noAGN) or a combination of galax- 
ies and AGN templates (AGN). 



6. Host galaxy properties 

6.1. Stellar mass 

In this section we estimate the stellar masses of the Type-2 
QSOs host galaxies. Several previous studies of Type-2 AGN 
host galaxy properties have assumed that the contamination of 
the AGN continuum to their measurements is negligible (e.g. 
'Silverman et al. "2008'; 'Silverman et alJ|2009l ISchawinski et al. I 
2010; Xue et al. 2010). However, those studies were focused on 
moderate luminosity AGNs, e.g. log(Lx < 44) erg s"^ . The AGN 
sample presented in this work consists of obscured nuclei, but 
they are all intrinsically powerful quasars, hence may still con- 
tribute in a non-negligible way to the overall photometry. It is 
then important to model this contamination, which would oth- 
erwise bias the measurement of the stellar mass from the SED 
fitting. To quantify the eff'ect of the nuclear component, we per- 
formed the SED fitting using either the combination of galaxy 
and AGN templates described in the previous section as well 
as only galaxy templates. In the top panel of Fig. O we com- 
pared the stellar mass estimates obtained with the two sets of 
templates. While for the majority of the objects, the introduction 
of an AGN component in the SED fit is irrelevant, there are a 
10% of cases where neglecting these components would cause 
more than a factor of two overestimation of the stellar masses. 

From the stellar masses reported in Table 1, one can im- 
mediately see that the hosts of Type-2 QSOs are mainly mas- 
sive galaxies. To quantify this impression, we need to compare 
the stellar mass distribution of these hosts with a representative 
sample of galaxies that matches their properties. We selected a 
parent sample of galaxies detected in the 3.6 jjm IRAC chan- 
nel from the COSMOS survey using the stellar masses mea- 
sured by Ilbe rt et aO (2010). We considered only the redshift 
range z=[0.8,1.5], where a high completeness in stellar mass 
above 5 x 10^ M© and accurate photometric redshifts are en- 
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Fig. 7. Stellar mass distribution of the Type-2 QSOs host galax- 
ies (solid histogram) and of the galaxies parent sample (dashed 
histogram) with 0.8 < z < 1.5. The second histogram has been 
rescaled for graphical reasons by a factor ~ 1000. The data 
points represent the fraction of galaxies hosting Type-2 QSOs 
in bin of stellar masses (see Sec. 16. II) . 



sured (see Fig. 8 of lllbert et al. II201Q|) . The final parent sample 
includes 76,061 galaxies. The sensitivity limits of the XMM- 
Newton observations vary spatially across the field (see Fig. 5 of 
ICappelluti et al. 2009), hence to derive the fraction of galaxies 
hosting a Type-2 QSO as a function of the stellar mass, we fol- 
lowed th e technique discussed in Sec. 3.1 of iLehmer et al. ll2Q07l 
(see also 'Silverman etalj|2008l |2009l). The Type-2 QSOs frac- 
tion for each stellar mass bin is determined by summing over the 
sample of Type-2 QSOs (Nqso2) 



A^QS02 



(1) 



where Ngai,i represents the number of galaxies for which we 
might have detected an AGN with X-ray luminosity Lx,i. We 
have plotted in Fig. [7] the fraction of galaxies hosting a Type-2 
QSOs . This fraction monotonically increases with stellar masses 
going from 0.02% at ^ 8 x 10^ M© to ^ 0.6% at ^ 10^^ 
Mo. 

This trend is consistent to what has been found for the over- 
all AGN population (Kauffmann et al. II2003I: iBest et al. II2005I: 



[Silverman et"Z1l2009TlBmsa et al. ll2009l:IXue et al. \\20m . 



For our sample of Type-2 QSOs, we do not detect broad lines 
and therefore we cannot derive "virial" black hole masses us- 
ing the cal ibration obtai ned with reverberation mapping on local 
AGN (e.g.l Kaspi et"aLl[2Q0Q) . To roughly estimate the mass of 
the black holes powering these obscured quasars, we used the lo - 
cal Mbh-M^ scaling relation provided bv Haring & Rix I (12004 . 
Since this relation was derived locally, we too k into account 
a_possible redshift evolution parametrized by iMerloni et al. I 
^10): Alog(MBn/M^)(z) = ^2/^g(l+z), where^2 = 0.68+0.12. 
We found that most (~ 90%) of the Type-2 QSOs have black 
hole masses in the range log(MBH) = [7.5, 9.5] M©. Using these 
estimates of black hole masses and the bolometric luminosities 
derived at the end of Sec. [21 we computed the Eddington ratio 
(^Edd = Lboi/^Edd, where LEdd = 1.38x10^^Mbh/Mo). We found 
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Fig. 8. Left panel. Color-mass diagram in the redshift range 
0.8 < z < 1.2 for the comparison sample of normal galaxies 
(gray contours), and the Type-2 QSOs hosts. The rest-frame col- 
ors of the Type-2 QSOs hosts have been measured after remov- 
ing the contribution from the AGN. Star-forming galaxies are 
indicated with the star symbol, while the circles correspond to 
passive galaxies. The dashed curve is the dividing line between 
re d passive and blue star-forming galaxies at z= 1 according 
to iPeng et aP (|201Q|) . The vertical bars show the variation in 
rest-frame colors if an unobscured AGN component is included 
in the SED fitting procedure (see Sec. 5). Central panel. U-B 
rest-frame color distribution for the Type-2 QSOs host galaxies 
(filled histogram) and the comparison sample of normal galax- 
ies (hatched histogram) in the redshift range 0.8 < z < 1.2. 
Right panel. As in the central panel but only for sources with 
> 3 X 10^0 Mq. 



that 50% of the Type-2 QSOs have /lEdd > 0.1, the remaining 
half being in the range 0.01 < /lEdd < 0.1. 



6.2. Rest-frame colors 

Galaxies display a co lour bimodali ty in b oth the local Universe 
(IStrateva et al. I boOll: iHogg et a l. ' "2002'; "Blanton et al. ' '2()03h 
Bell et a l. 2004, Weiner et al. 2005rup 



and at hi gher redshift 
to 1 : iFranzetti et al 



2007, Cirasuolo et al. 2007 up t o z^ 



1.5:iGiallongo et al. 2005. Cassata et al. 2008, WilHams et al. I 
Hoof up to 2). This bimodality in red-sequence and blue- 
cloud galaxies has been interpreted as evidence of a dichotomy 
in the star-format ion and merging histories of galaxies (e.g. 
iMenci et al. Il2005l) . Color-magnitude diagrams have been used 
as tools in galaxy evolution studies, and since many models in- 
voke AGN feedback as an important player in this evolution, it 
may be interesting to locate the hosts of Type-2 quasars in those 
diagrams. 

From the galaxy component of the best fit to the Type-2 QSO 
SED (see Sec. O, we derived rest- frame colors for the hosts 
in the Johnson-Kron-Cousins system. The inclusion of an AGN 
component in our SED fitting procedure prevents contamination 
from the nuclear emission to the colors of the host galax}|3. 

From the histogram in the central panel of Fig. [8] 
one could conclude that the host galaxies of Type-2 
QSOs mainly have rest-frame colors in-between the blue 
and red galaxy populations, similarly to what has al- 
ready been reported for lower luminosity X-ray /optically 



^ See the discussion in Appendix A of lHickox et al. I (l2009h about the 
contribution of the nucleus to optical colors of X-ray and radio- selected 
AGNs if the SED fitting is done using only galaxy templates. 
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selected AGNs (iNandra et al. I I2QQ7I: Martin et al. I l2007t 
'Rovilos & GeorgantouDoulos " "2007"; 'Silverman et al. " "20081 
Georgakakis et al. 2008; Schawinski et al. 2009; Hickox et al!] 
120091) . ISilverman et al. I (120091) demonstrated that if one consid- 
ers stellar-mass rather than luminosity selected samples of AGN 
and parent galaxies, the fraction of ga laxies host i ng an AGN 
does not decline towards bluer colors. IXue et al. I (l2010l) con- 
firmed that the AGN fraction remains nearly constant as a func- 
tion of color for a stellar mass selected sample of AGN from the 
Chandra Deep Fields. Both these works focused on AGN with 
luminosities below Lx < 10"^"^ erg s"^. Our sample of Type-2 
QSOs is by construction brighter and therefore care is required 
when extrapolating these previous results. If we select only the 
objects with > 4 x lO^^M© f or which our COSMOS sam- 
ple is mostly complete (see Il bert et al. ^'2010^. the distribution 
of rest-frame colors is modified as shown in the right panel of 
Fig. [HI In a mass-selected sample of galaxies the distribution of 
rest- frame colors for active and inactive galaxies is similar. Since 
stellar mass is likely playing a major role in determining the evo- 
lution of galaxies, mass-selected samples are more suitable when 
comparing the properties of the AGN hosts with those of inactive 
galaxies. 

If we adopt the dividing line bet ween blue star- form ing and 
red passive galaxies as suggested by lPeng et al. I (1201 Ol) 



(U-B)ab = IAO + O.015log( 



lOiOMc 



-)-0.18z. 



using z = 1.0, we find that ^ 28% of the host galaxies of 
Type-2 QSOs in the redshift range 0.8 < z < 1.2 are red and 
^ 7 2% are blue. If \ yere to use the (U-B) versus Mb diagram and 
the IWillmer et al. I (l2006l) color division, the fraction of Type- 
2 QSOs hosts that are classified as red galaxies would increase 
to ^ 58%. We stress that the use of rest- frame colors to study 
the properties of AGN host galaxies is subject to a proper treat- 
ment of reddening. Objects with intermediate or red colors could 
actually be star-forming g alaxies obscured by large amounts of 
dust on kpc scales (e.g. see lBrammer et al. ]|2003 Pozzetti et al. 
120101) . As we describe in the next section, we were able to de- 
rive estimates of the SFR of our host galaxies from the SED 
fitting. We can therefore use the specific star-formation rate 
(SSFR=SFR/M^) to divide the sample into active and quiescent 
galaxies depending on whether log(SSFR/Gyr"^) is above or be- 
low -1, i.e., galaxies that would take less or more than 10 Gyr 
to double their stella r mass at the present SFR, respectively (see 
iPozzettietaD 120101) . We find that ^ 20% of the red galaxies 
are actively forming stars, and therefore their colors are due to 
the presence of dust rather than an old stellar populatio n. Our 
result is in agreement with previous studies. Cardamone et al. I 
(I2OIOI) . using rest- frame near-infrared colors of lower luminos- 
ity AGN hosts to separate passive from dusty galaxies, found 
that ^ 2 5% of the hosts in th e red cloud are dusty star-forming 
galaxies: ferusa et al. I (l200 9b) from the SED fitting of a sample 
of obscured AGN in the Chandra Deep Field South reached the 
conclusion that the fraction of dusty star-forming hosts was up 
to ^ 50% of their sample. 



6.3. SFR 

As described in Sec.O we estimated the SFR values of the host 
galaxies using the SED fitting technique. We stress that to reli- 
ably estimate the SFR of the host galaxy it is crucial to remove 
the contribution of the AGN to the observed photometry. As al- 
ready discussed, we achieved this by performing the SED fitting 
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Fig. 9. SFR versus stellar masses obtained from the SED fitting 
for the host galaxies of the Type-2 QSOs in the redshift range 
0.8 <z< 1.2. Galaxies with log(SSFR/Gyr-^)> -1 are indicated 
with a star symbol, otherwise a circle was used. Type-2 QSO 
hosts with red (U-B) rest- frame colors according to the color 
division by IWillmer et al. I (l2006l) have been highlighted with a 
smaller lighter circle. The vertical bars show the variation in SFR 
if an unobscured AGN component is included in the SED fit- 
ting procedure (see Sec. 5). The dashed lines are the correlation 
and its one-cr confidence intervals found for blue star-forming 
galaxies by H lbaz et al.J (l2007l) and PACS/Herschel sources by 
Rodighiero et aLj (i2010 ) at z^^ 1 . The curves have been rescaled 
from the original Salpeter IMF to the Chabrier IMF adopted in 
this work. The vertical dashed line indicates the mass complete- 
ness limit of the PACS sample. The triangle indicates the av- 
erage SFR for the hosts of X-ray selected AGNs estimated by 
iLutz et al. I (l2010h by stacking 870 jim LABOCA/APEX data. 



with a combination of AGN and galaxy templates. To quantify 
the effect of not modeling the nuclear emission, we compared 
the SFR estimates obtained only with galaxy templates and a 
combination of galaxies and AGN SEDs, as discussed in Sec. 
[6Jl We find that for ^ 40% of the objects neglecting the AGN 
component will overestimate the SFR by more than a factor of 
two, and up to a factor of fifty for a few galaxies (see bottom 
panel of Fig. [6]). 

For comparison, we also estimated the SFR us ing the UV 
continuum luminosity. Using the prescriptions of iDaddi et al. I 
(120041) for star-forming galaxies at 1.4 <z< 2.5 and assuming a 
Chabrier IMF, we have 



jcorr _ J ^ y in0.4xAi5oo 
^1500A - ^1500A X 



SFR(MQyr-^) = 



where L^5qq^ is the galaxy luminosity at 1500 A (erg/s/Hz), 
MsooA ^^^^^^ted for extinction by dust, and A 1500 ~ 10 xE(B- 
V) assuming a Calzetti's law for the dust extinction, where we 
use the E(B-V) values derived from the best SED fit using both 
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Fig. 10. Evolution of the average SSFR for Type-2 QSOs hosts 
as a function of redshift. For each redshift bin we report the aver- 
age stellar mass of the host galaxies. Empty circles indicate the 
average SSFR per redshift bin if an unobscured AGN component 
is included in the SED fitting procedure (see Sec. 5). The contin- 
uum line reproduces the evolution law found by iPannella et al. I 
(i2009i) for star-forming galaxies up to z~ 2. The dashed line in- 
stead is t he b est fi t to the SSFR redshift evolution derived by 
iKarim et al. l (l2Qllh in the mass range log(M^)=[10.6,11.0] M©. 
The two squ ares are the SSFR for n ormal galaxies at z= [4, 5] as 
presen ted by lGonzalez et al. I (12010) using data from lStark et al. I 
(l2009h 
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Fig. 11. Evolution of the average SSFR for Type-2 QSOs hosts 
with log(M^)=[10.6,11.0] Mo from z~ 1 to z~ 2. The lines and 
symbols are the same as in Fig.fTOl 



galaxies and AGN templates. Con sidering the ^ 60 QSO-2 hosts 
that match the lDaddi et al. I (12004 ) selection criteria, we find that 
the median ratio (SFR5£/)/SFRf/y) is ^ 1.1 with a dispersion of 
^0.6. 

One of the main goals of this work was to compare the level 
of star formation in the hosts of the Type-2 QSOs with that mea- 
sured for ''non-acti ve" galaxies. Surveys of the local Universe 
(Brinch mann et al. 11200 4: Salim et al. 2005) found that galaxies 
can be divided mainly into two populations: star-forming objects 
exhibiting a continuous increase in SFR with along a tight 
correlation and quenched galaxies with little star formation and 
high stellar masses. The tight correlation between SFR and 
for star-forming galaxies has been confirmed to a lready be in 
place at higher re dshifts (iNoeske et al. ll2QQ7a[ andlElbaz et al. 
120071 up to 1; iDaddi et al. I I2QQ7L and IPannella et al. 1 12009 
up to z^ 2), also whe n using PACS/Herschel observations 
(iRodighiero et al. IbOlQ ). We refer to it as the main sequenc e 
(MS) of star-forming galaxies following iNoeske et al. I (l2QQ7ah . 
The slopes of the MS at diff'erent redshifts appear to be the same 
but the normalization has decreased by a factor ^3.7 going from 
z^ 2 to z^^ 1 , and an additional factor 7 at z^^ ( Daddi et al. I 
I2QQ7) . This decrease in the activity of star-forming galaxies of 
a given stellar mass reflects the decline in the global SFR den- 
sity with time, and a likely cause of this is gas exhaustion (e.g. 
"Noeske et al. "2007b'). 

The galaxy dichotomy found using the correlation between SFR 
and M* is more fundamental than that based on rest- frame col- 
ors: red colors do not necessarily mean that the galaxy has 
stopped forming stars (see the discussion in Sec. l6.2| and the red 
colors of some MS galaxies at high M* in Fig. 17 of El baz et al. I 
120071 Our aim is to compare the properties of the host galaxies 
with a parent sample of normal galaxies in the context of this 
correlation. 

In the redshift range 0.8 <z< 1.2, we have 35 Type-2 QSOs; 
and we were able to establish a good fit to their SED for 340 
of which. We find that ^ 62% of the Type-2 QSO host in this 
redshift range are actively forming sta r s. In t he sample of star- 
forming galaxies used by ' Elbaz et aTl (l2QQ7l) to derive the MS 
dii 1, only "blue" galaxies were considered (see their Fig. 
13). For consistency, we derived rest- frame absolute magnitudes 
in the standard U and B Johnson's filters from the galaxy com- 
ponent of the SED fit of our Type-2 QSO hosts, and separated 
them into blu e and red galaxies using the color division by 
IWillmer et al. I ( |2QQ6| ). As shown in Fig. |9j there is a mild cor- 
relation between SFR and M^ for the "blue" star-forming hosts. 
The Spearman rank correlation coefficient for these sources is 
r= 0.38 (p= 0.2). For the majority of our "blue" star-forming 
host galaxies, the level of SSFR at 1 is consistent with the 
MS of star-forming galaxies at the same red shift. In Fig . |9l we 
also indicated, the average SFR estimated by lLutz et al. I (l201Qh 
for the hosts of X-ray selected AGN at z~ 1 obtained by stacking 
870 yum LAB OCA observations in the Extended Chandra Deep 
Field South. Most of these hosts have stellar masses in the range 
10.4 <log(M^) < 11 Mo (Xueetal. 2010). From this stacking 
analysis, the derived value of SSFR for AGN hosts is also con- 
sistent with the MS of star-forming galaxies at z~ 1. 

We repeated a similar analysis for the hosts in the three 
higher redshifts bins: 1.2 < z < 1.8 (36 objects), 1.8 < z < 2.5 
(28 objects), and 2.5 < z < 3.2 (12 objects). First, we find that 
the fraction of star-forming galaxies between our Type-2 QSO 
hosts increases with redshift: 66%, 71%, and 100% respectively. 



^ Source 5326 has near-lR photometry contaminated by a nearby 
source and will not be considered in this analysis. 
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Fig. 12. Type-2 QSOs hosts classified as bulge-dominated gaisix- 
ies. The images are 10'' on a side. They were obtained by com- 
bining Subaru B-band (blue), ACS F814W (green), and CFHT 
Ks (red) images. XID and coordinates of the source are reported 
at the top of each image. 



We then looked at the evolution of the SSFR for Type-2 QSO 
hosts between redshift z= 0.8 and z= 3 as shown in Fig.[TOl At 
a first look, this figure would suggest that there is only a mild 
evolution, if any, of the average SSFR for Type-2 QSO hosts go- 
ing from z~ 1 to z~ 2, while there is no evolution at all at z> 2. 
The latter agrees with the constant SSFR measured from z~ 7 
to z~ 2 for normal galaxies (e.g. Gonzalez et al. 2010), while 
the former disagrees wi th the substantial dro ps o f the average 
SSFR below z~ 2 (e.g. iNoeske et al. ll2007aL and lDaddi et"^ 
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Fig. 13. Type-2 QSOs hosts classified as disk-dominated galax- 
ies. The images are 10'' on a side. They were obtained by com- 
bining Subaru B-band (blue), ACS F814W (green), and CFHT 
Ks (red) images. XID and coordinates of the source are reported 
at the top of each image. 
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Fig. 14. Type-2 QSOs hosts classified as possible-merger. The 
images are 10" on a side. They were obtained by combining 
Subaru B-band (blue), ACS F814W (green), and CFHT Ks (red) 
images. XID and coordinates of the source are reported at the 
top of each image. 



I2007h . IShao et al. I (l2010l) used Herschel/R\CS observations of 
the GOODS-N field to estimate the star-formation rate in X-ray 
selected AGNs. They find that the host far-infrared luminosity 
of AGN with L2-meV ~ 10"^^ erg s"^ does not significantly in- 
crease from z~ 1 to z~ 2 (see their Fig. 4). Since the rest- frame 
far-infrared emissi on is typically dom inated by the host rather 
than the AGN (e.g. Netzer et al. "2007b) it can be used as a SFR 
diagnostic and therefore the Shao et al. (2010) results are con- 
sistent with ours in Fig.[TOl However, submillimetre surveys of 
narrow line radio galaxies, possibily radio loud Type-2 QSOs, 
found that their FIR luminosity rises with redshift out to z^ 4, 
and interpreted it as an increase in the star formation activity 
with redshift (e.g. Archibald et al. 2001; Reuland et al. 2003). 
We note that the samples of Type-2 QSOs in the redshift bins 
plotted in Fig. [TOl have diff'erent mean stellar masses: at z~ 1, 
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Fig. 15. Location of the hosts morphologically classified as 
bulge-dominated (filled circles), disk-dominated (filled ellipse), 
possible-mergers (filled stars) in three planes defined using 
the four non-parametric measurements of the light distribution: 
asymmetry, conce ntration, M2Q, Gini c oefficient. The contours 
include 99% of the lScarlata et al. I (I2QQ7I) galaxies with morpho- 
logical classification as bulge-dominated (red), disk-dominated 
(blue), and possible-mergers (green). 



the average stellar mass is < >~ 3 x 10^^ M©, while it 
is high er, < >~ 6 x 10^^ M©, at z~ 3. Also the curve 
from iPannella et al. I (I2QQ9 ') plotted in Fig. [TOlis niost va lid for 
~ 3 X 10^^ Mo. According to iKarim et al.] (|201l|) , there 
is a correlation between SSFR and stellar masses, with higher 
mass objects having lower SSFRs. For the three redshifts bins 
below z~ 2, we now include only the hosts with log(M-A.)=[10.6- 
11.0] M© and consider the best fit to the redshift evolution of 
th e SSFR for norma l galaxies in the same mass range (Table 5 
of lKarim et al. 1120111) . As Fig.Oshows, the agreement between 
Type-2 QSO hosts and normal galaxies is now excellent: the av- 
erage SSFR for our hosts decreases by a factor 3.6 going from 
z~ 2 to z~ 1. 



7, Morphology 

Another important diagnostic of galaxy evolution is their mor- 
phologies, which correlates with their star-formation histories. 
In Type-2 QSOs, obscuration on parsec scales prevents the light 
of the nuclei from outshining the host galaxy, therefore mak- 
ing it easier, compared to unobscured Type-1 QSOs, to study 
the morphological properties of their hosts. Using the Cycle 12 
Hubble Space Telescope (HST) Advanced Came ra for Surveys 
(ACS ) F814W images of th e COSMOS field ( Scoville etal] 
12007 1: iKoekemoer et aDl2007l), we studied the morphology of 34 
Type-2 QSO host galaxies with z < 1.2 and Tab < 24. We used 
an upgrade d version of the Z urich Estimator of Structural Types 
(ZEST; Sca rlata et al.1l2007l) known as ZEST-F (Carollo et al, 
in preparation). This new tool, based on a principal component 
analysis, includes additional measurements of non-parametric 
morphological indices for characterizing both structure and sub- 
structures. For each of the 34 host galaxies, we measured the 
half-light radius (ri/2), the Petrosian radius (Rp), and the el- 
lipticity (5). In addition, we derived estimates of four popular 
non-par ametric measurements of the light distrib ution in galax- 
ies (e.g. [Abraham et al. Il2003t iLotz et al. Il2004 : concentration 
(C), asymmetry (A), second-order moment of the brightest 20% 
of the galaxy flux (M20), and Gini coefficient (G). The values of 
these parameters for each of the 34 Type-2 QSOs hosts are re- 
ported in Table 2. Using a combination of this information and 
after a visual inspection of the ACS images, we divided our host 
galaxies into three classes: bulge-dominated, disk-dominated, 
and possible-merger. The breakdown in morphological classes 
for the Type-2 QSO host galaxies at z< 1.2 is the following: 
57% bulge-dominated, 23% disk-dominated, and 20% possible 



mergers. Color cutouts of the hosts divided by morphological 
classes are shown in Fig. [121 El and [141 In Fig. [151 we plot 
the distribution of our host gala xies in three planes: G vs. M20, 
C vs. M20, and C vs. A. IScarlata et al. I (l2007l) obtained reUable 
morphological classifications for ~ 56000 Iab < 24 galaxies 
in the COSMOS field using ZEST. For comparis on, we added 
in Fig. [151 the contours enclosing 99% of the Scarl ata et al. I 
(2007) galaxies belonging to the ZEST morphological classes 
T=l (red), T=2.0 or 2.1 (blue), which correspond to disk galax- 
ies with a Sersic index n > 1.25, and T=3 (green). While in 
the first two panels our morphological classificatio n agrees quite 
well with the large galaxies sample studied by Scarlata et al. I 
(2007), in the third panel a significant fraction of our hosts 
are outside the contours mainly due to larger asymmetry val- 
ues than observed for galaxies of the same morphological class. 
This could be interpreted as an indication of a higher fraction 
of disturbed morphologies in our sample of Type-2 QSO host 
galaxies than normally observed in a magnitude limited sam- 
ple of normal galaxies. Furthermore, our hosts tend to be highly 
concentrated and ^ 65% of them have C> 3, a range in which 
elliptical and bulge-dominated galaxies are located. This is not 
surprising since they are mostly massiv e galaxies that usually 
have a consistent bulge already in place. S ilverman et al. I (l2008k 
also found that for X-ray selected AGN of lower luminosity (log 
Lo.5-iokeV < 43.7 erg s"^) at z< 1 the majority of the host galax- 
ies are bulge dominated. We note that in the morphological anal- 
ysis presented in this paper, we did not attempt to subtract any 
possible residual point-like source at the center of the galaxy. If 
present, these point sources wou ld bias the co ncentration mea- 
surements towards higher values: iGabor et aTl (1200 9) found that 
this is the case for a sample of ^ 200 X-ray selected AGN at 
0.3 <z< 1 in the COSMOS fi eld (see their Fig. 12). The AGN 
sample of iGabor et al. I (l2009l) was, on average, less obscured 
than the Nh > 10^^ cm"^ objects we present, hence the like- 
lihood that a residual point-like source is present in the ACS 
images is higher. We plan to study this possible bias in detail in 
future work. 

As reported above, we were able to detect signs of recent merg- 
ers in only ^ 20% of the hosts of Type-2 QSOs at z< 1.2. The 
lack of morphological merger signatures for the majority of the 
objects should be treated with caution. The ACS images could 
simply be not deep enough to detect faint surface brightness fea- 
tures, as tidal tails or shells, to trac e past interactions. For ex- 
amplel^mos Almeida et al.l (|201l|) . using deep optical images, 
found that ^ 80% of a complete sample of z< 0.7 narrow line 
radio galaxies show signs of interactions based on the detection 
of faint morphological features. In addition, some models (e.g. 
"Pi Matteo et al. I l2(}05b predict that the quasar phase occurs to- 
wards the end of the merging process, and at that point the final 
galaxy sho uld be most ly relaxed and have low values of asym- 
metry (e.g. IConselice1 f2003). In this scenario, it becomes diffi- 
cult to test merging-induced quasar fueling models by studying 
the morphological appearance of the host galaxies. At the same 
time, this hypothesis also opens another possible approach: a 
major merger would most likely form spheroidal galaxies and 
large disks should usually be destroyed 0. As described above, 
we find that the majority of the Type-2 QSO hosts at z< 1.2 



^ Some s imulations (e.g. iBames & HemquisTI 119961 : 
ISpringel & He mquist 2005) show that mergers of gas-rich galax- 
ies could produce a disk-dominated remnant but only for particular 
orientations of the angular momentum vectors. These conditions are 
usually not related to substantial BH growth (iHopkins & Hemquist I 
i200a). 
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Fig. 16. Distribution of Eddington ratios for the hosts of Type- 
2 QSOs morphologically classified as bulge-dominated (top 
panel), disk-dominated (middle-panel), and possible mergers 
(bottom panel). 



are bulge-dominated galaxies and that a few of them contain 
disks. In a study of the morphology o f the hosts of 14 AGN 
in the XMM-COSMOS survey, Ciste mas et aTl (l201ll) found 
that: a) there is no significant diff'erence in the distortion frac- 
tions between active and inactive galaxies; b) ~ 65% of the 
AGN hosts are disk dominated. On the basis of these findings, 
they concluded that major mergers are not the trigger of B H ac- 
cretion . We note that the vast majority of the Ciste mas et al. I 
(I2OIII) AGN sample has lower X-ray luminosities than the Type- 
2 QSOs presented in this paper: only 18 sources (13%) of their 
sample have Lx > 10"^^ erg s"^ and these are mostly unobscured 
QSOs. The diff'erent fraction of disk galaxies between the two 
studies could be connected to diff'erent triggering mechanisms 
of accretion activity at lo\y and high AGN luminosities (e.g. 
[Hopkins & Hemquist""2009'). 

In Fig. [161 we compared the distribution of Eddington ratios 
for the three morphological classes. The range of values of this 
statistic is limited, but bulge-dominated galaxies tend to host 
AGN with low Eddington ratios (A < 0.1), while disk-dominated 
or merging galaxies at their centers contain BHs accreting at 
high Eddington ratios (A > 0.1). The mean BH mass of the 
bulge-dominated sample is < Mbh >~ 4 x 10^ M©, while for 
the mergers < Mbh >~ 1 x 10^ M©. The distributions in Fig. 
[16] are therefore consistent with the ge neral finding that the low- 
est m ass BHs are the fastest accretors (IMcLure & Dunlopl2002l: 
iN^tzer & Trakhtenbrot.2007.) . 

8. Conclusions 

We have explored the relationships between SMBH accretion 
in obscured quasars and the physical properties of their host 
galaxies. The sample of 142 Type-2 QSOs was selected from 
those considered for the X-ray spectral analysis in the XMM- 
COSMOS survey. The bulk of this quasar sample is located in 



the redshift range 1 < z < 2.5, which includes the peaks of both 
the global star-formation rate and AGN activity. The obscured 
quasars presented in the paper have a mean bolometric luminos- 
ity of Lboi = 8 X 10^^ erg s"^ a regime where several models 
invoke major mergers of gas-rich galaxies as the main fueling 
mechanism of the central SMBH. 

A key ingredient of this work has been the ability to use an 
AGN-hhost galaxy SED decomposition technique to put reliable 
constraints on the stellar masses and star-formation rates of the 
hosts. This has been possible for the excellent multi- wavelength 
coverage in the COSMOS field. 
The main results of this work are the following: 

- Type-2 QSOs reside almost exclusively in massive galax- 
ies: more than 80% of their hosts have M* > 10^^ M©. 
The fraction of galaxies hosting Type-2 QSOs monoton- 
ically increases with stellar masses going from 0.02% at 
M^ 8 X 10^ M© to ^ 0.6% at M^ ^ 10^^ M©. 

- The rest-frame colors of their host galaxies are mostly be- 
tween those of the blue and red galaxy population. We con- 
firm that this is the eff'ect of a luminosity-selected rather than 
a mass-selected sample. Furthermore, the classification of 
the host galaxies based on their rest- frame colors is subject to 
the eff'ect of dust reddening: ^ 20% of the host galaxies with 
red or intermediate rest- frame colors are dusty star-forming 
objects. 

- At z~ 1, ^ 62 ± 7% of Type-2 QSOs are actively forming 
stars and this fraction increases with redshift, 71 ± 3% at 
z~ 2, and 100% at z~ 3. 

- We find that the redshift evolution of the specific star forma- 
tion rate for the hosts of Type-2 QSOs is in excellent agree- 
ment with that observed for star-forming galaxies in the red- 
shift range 1 < z < 3: the average SSFR increase by a factor 

3.6 going from z~ 1 to z~ 2, and remains constant at 
higher redshifts. 

- A morphological analysis of the ACS images has indicated 
that the majority of the host galaxies are bulge dominated, 
there being a few cases of disk galaxies or mergers. This 
lack of merger signatures cannot be automatically translated 
as proof that mergers are not playing an important role in 
the accretion history of these quasars. Several models in- 
deed assume that the quasar phase occurs toward the end of 
the merging process, and at that point we expect the host 
galaxy to be mostly relaxed with low values of asymme- 
try. Numerical simulations also predict that large disks are 
destroyed during a major merger, therefore detecting large 
number of disks could place some strong constraint on the 
occurrence of a major merger. With the current data, we have 
detected clear disks in only a few cases. 

- We find that bulge-dominated galaxies tend to host Type- 
2 QSOs with low Eddington ratios (A < 0.1), while disk- 
dominated or merging galaxies have at their centers BHs ac- 
creting at high Eddington ratios (A > 0.1). 

The morphological analysis presented in this paper is based 
on F814W ACS images, which corresponds to the rest- frame 
UV at the average redshift of our galaxies. This implies that we 
are mainly mapping regions of recent star-formation. It will be 
extremely interesting to repeat this analysis using future near- 
IR observations from space that will provide high resolution 
images in the optical rest-fram e. The CANDELS HST Multi- 
cycle Treasury Program (Grogin et al. [[201 1 ; Koekemoer et al. L 
1^1 ) will image with WFC3/HST part of the COSMOS field 
and therefore provide a less biased view of the morphology of 
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these Type-2 QSOs hosts. It would be also interesting to es- 
timate the nebular continuum contribution to the UV SED of 
Type-2 quasars host galaxies. This would be possible using NIR 
spectroscopy to derive Hyg flux and the line of sight reddening 
for each quasar (see e.g. lTadhunter et aini2002h . Finally, the on- 
going Herschel survey of the COSMOS field, as well as other 
patches of the sky, will improve the mid and far-infrared cov- 
erage, a crucial wavelength range to measure the SFR of AGN 
hosts. 
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Table 1. SED of Type-2 QSOs. The legend for the column named Survey is the following: s= spectroscopic redshift; p=photometric redshift. 
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The photometry is contaminated by the halo of a bright nearby star 
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Table 1. continued. 
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L0.8 
L0.8 


/TC c65.5 


OA "2 

-20.3 


01 A 

-Z1.4 




AA 

44. JU 


00 n23.4 
23.O225 


1 onn 


p 


11.0 


Ll.O 
Ll.O 


^•^0.8 


on Q 


00 1 
-ZZ. 1 


5054 


44.54 


22.022-3 


1.420 


p 


10.8 


L0.8 
L0.8 


38-9^^:^ 


-19.8 


-21.0 


j0j6 


A A 1 O 

44. ly 


00 c22.8 
^^"^22.3 


1 000 
l.ZZU 


p 


11.0 


Ll.O 
Ll.O 


n 70.7 


OA n 

-20. / 


010 
-21.9 


5057 


44.19 


22 2 
"^"^•"^20.8 


1.690 


p 


11.1 


Ll.l 
Ll.O 


1 1 r" 1 1 1 6 7 

115 1 

^^^•^115.1 


-21.5 


-22.5 


5066 


44.42 


22 422-'7 

z.z..^2i.9 


0.970 


p 


11.0 


Ll.O 
Ll.O 


9 5^-6 


-20.8 


-22.0 


5070^^ 


44.77 


00 o23.4 
z.j».z.23j 


1.450 


p 


10.3 


L0.3 
L0.3 


462 4"^^^-^ 

^^^^.^461.3 


-24.7 


-24.9 


5111 


44.63 


22.9i:^ 


1.080 


p 


10.0 


LO.l 
LO.O 


12.411: 


-19.7 


-20.6 


5120 


44.98 


22 422-6 


3.650 


p 








0.0 


0.0 


5149 


44.07 


00 o22.8 
^^•"^20.4 


1.560 


p 


10.6 


10.7 
10.6 


93-8^^:^ 


-20.6 


-21.5 


5153 


44.13 


22 923-1 


0.787 


s 


10.0 


LO.O 
0.0 


47 1^^-' 


-19.6 


-20.5 



The ACS image shows a source at 0.5" from the chosen counterpart, therefore the IR photometry is the blending of the emission 
sources. 

For this objects we have photometry only in five filters. 
1^ The photometry is contaminated by a nearby bright star. 
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Table 1. continued. 



XID 


log Lx[0.5-10keV] 


log Nh 


z 


Survey 


log M* 


SFR 


Mu 


Mb 




(ergs-^) 


(cm-2) 






(Mo) 


(Mo yr-i) 


(AB) 


(AB) 


5158 


44.22 


99 422.6 
^^•^22.0 


1.440 


P 


10 9^0.9 
^'^•^10.9 


22.9i:^ 


-21.7 


-22.5 


5161 


45.54 


22 Q^^-^ 


3.160 


P 


11 1^1-^ 


177.4}77-4 


-22.6 


-23.4 


5162 


44.94 


22 522-9 


3.524 


s 


10 5^0.5 


699 8^99.8 

'^^^•^695.0 


-23.4 


-23.8 


5178 


44.43 




1.720 


P 


11 3^1-^ 


1 31-^ 


-21.4 


-22.7 


5186 


45.16 


22 622-7 


2.540 


P 


Q q9.9 
^•^9.9 


195.4^gg g 


-20.8 


-21.4 


5229 


44.37 


22 623-1 


1.030 


P 


10 7^0.8 


11 III-7 


-19.8 


-20.9 


5232 


44.73 


22 923-2 


2.830 


P 


11-211:? 


79 o89.7 


-20.6 


-21.7 


5266 


44.72 


23 123-3 
^^•-^22.8 


1.950 


P 


10-81^:^ 




-20.7 


-21.7 


5278 


44.94 




1.620 


P 


10 6^0-^ 


32.61:^ 


-20.7 


-21.6 


5288 


44.15 


99 922.4 


1.000 


s 


11 lii-i 


2 22-2 


-20.7 


-21.9 


5292 


44.57 


99 423.0 

"^"^•^20.4 


2.730 


p 


11 O^^-O 


27.0^^:^ 


-19.7 


-21.0 


5305^ 


44.17 


99 423.0 

"^"^•^21.4 


2.500 


p 




161-81^^:^ 


-21.4 


-21.9 


53262i 


44.35 


22 723-2 
^22.0 


0.950 


p 




9 29.3 


-19.3 


-19.6 


5347 


44.65 


23 124-3 
"^^•^20.4 


3.120 


p 


11 411-4 


212.8^4.9 


-22.2 


-23.2 


5357 


44.83 


23 623-9 


1.670 


p 


10 710-7 


1-61:^ 


-20.4 


-21.4 


5368 


44.17 


99 922.4 


1.230 


p 


11 2II-2 


105 9105-9 

^'^^•^105.9 


-20.9 


-22.1 


5369 


44.16 


22.1^?:^ 


0.950 


p 


10.4}0-4 


21 121-7 
"^^•^21.0 


-19.7 


-20.7 


5400 


44.19 


22.7i:^ 


0.680 


p 




1 2I-2 


-20.0 


-21.2 


5414 


44.10 


22 422-5 

"^"^•^22.2 


0.910 


p 


10 710.V 


O9I-O 


-19.7 


-20.9 


5417 


44.62 


99 422.8 


2.050 


p 


11 0^1-0 


115 3 51153-5 

^^^^•^1153.5 


-24.1 


-24.6 


5427 


44.86 


23 623-8 


1.177 


s 


11 lll-l 


272-7 


-21.5 


-22.6 


5441 


44.11 


22.1^?:^ 


1.170 


p 


11 1^1-1 


4 2^-2 


-20.8 


-21.9 


5448 


45.30 


22.81:° 


2.060 


p 


11 0^1-0 


43 0^3-1 


-21.8 


-22.7 


5452 


44.40 


22.022.4 


1.920 


p 


10 0^0.0 


131 51^^-0 

^^^•^131.5 


-21.2 


-21.7 


5467 


44.14 


90 924.6 

^■^•^20.4 


0.870 


p 


10 7^0.7 


1-61:^ 


-19.1 


-20.3 


5488 


44.36 


23 123-6 
^'22.2 


1.280 


p 


11 0^1-0 


163.7lS^ 


-21.8 


-22.8 


5496 


44.60 


23.6|:0 


0.694 


s 


10 610-6 
^'^•"10.6 


54 554-5 


-20.6 


-21.6 


5512 


44.31 




1.130 


p 


11 411-4 


5 2^-2 


-22.2 


-23.3 


5546 


44.62 


22.8204 


1.200 


p 


11 lii-i 


242.4 


-21.3 


-22.5 


5547 


44.07 


99 422.6 
"^"^•^22.0 


0.960 


p 


9 59-5 


17 517-8 


-19.2 


-19.9 


5548 


45.01 


23 123-5 

"^^•-^22.8 


1.060 


p 


10 910-9 


1 Oi-o 


-20.3 


-21.5 


5553 


44.51 


23 1-1-0 

"^■^•^20.4 


2.390 


p 


11 2II-2 


1 111 


-20.6 


-22.0 


5563 


44.88 


22.8i:« 


3.410 


p 


10 310-^ 

^'^•■^10.3 


47Q 9493.2 
■^•^443.6 


-21.7 


-22.3 


5569 


45.04 


22.7i:2 


1.810 


p 


11 2II-2 


3.3^-^ 


-21.1 


-22.3 


5573 


44.88 


99 922.4 
"^"^•"^22.0 


1.660 


p 


11 2II-2 


182.0}^^:^ 


-21.4 


-22.5 


5582 


44.60 


22 923-I 


1.490 


p 


10 710-7 




-20.3 


-21.4 


5606 


44.73 


22 623-2 

20.8 


4.166 


s 


10 210-2 

10.2 


401.8401-B 


-23.4 


-23.7 


10094 


44.69 


22 923-2 

"^"^•^20.4 


2.670 


p 


10 710-7 


5 0^-0 


-20.7 


-21.8 


10499 


44.11 


22 523-0 


1.990 


p 


11 2II-2 


98.99B-9 


-20.9 


-22.0 


10690 


44.52 


22 723-3 


3.100 


p 


10 210-2 

^^•^10.2 


177.8}79-i 


-22.7 


-23.1 


10719 


44.68 


22 522-9 
"^^•^22.1 


1.700 


p 


10 510-5 


360.63^4.8 


-22.9 


-23.3 


30292 


44.37 


23 0-^-0 


2.550 


p 


10 Olo-o 


1.1}-} 


-19.7 


-20.6 



The quality of the SED fit is poor. 

The IR photometry is contaminated by a nearby source. 
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Table 1. continued. 
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XID 


logLx[0.5-10keV] 


log Nh 


Z 


Survey 


log M* 


SFR 


Mu 


Mb 




(ergs-^) 


(cm-2) 






(Mo) 


(Mo yr-i) 


(AB) 


(AB) 


30361 


44.19 


22.8^:^ 


1.910 


P 


11 l^ii 


156 O'^'^o 


-21.2 


-22.2 


30789 


44.99 


22.8^3-3 


1.810 


P 


11 liii 


^•'2.7 


-20.9 


-22.1 


31357 


44.63 


23.222*9 


1.700 


P 


11 411-4 


1.71:^ 


-21.7 


-23.0 


31419 


44.35 


23 6-1-^ 


1.130 


P 


Q QlO.O 
^•^9.9 


14.6iy 


-18.6 


-19.6 


53351 


45.80 


23 2"^-^ 


3.000 


P 






-21.1 


-21.9 


54316 


44.26 


22 5^2.9 

"^^•^22.0 


1.980 


P 


11 211-2 




-21.1 


-22.2 


54514 


44.44 


23 6-1-^ 
^■^•"20.4 


0.707 


s 


11.111:1 


252-5 


-20.8 


-22.0 


60075 


44.15 


^^•^22.4 


1.080 


P 


10.8l°i 


1 3I-3 


-19.4 


-20.7 


60144 


44.07 


23.0->:« 


1.500 


P 


10 91^-^ 


5^-5 
^•^0.5 


-21.1 


-22.2 


60168 


44.84 


22.6^:^ 


2.140 


P 


10 91^-^ 


380.23^3.v 


-21.9 


-22.7 


60214 


44.17 


22.8-10 


1.640 


P 


10 91^-^ 




-19.8 


-21.1 


60255 


44.15 


22.8-o';0 


0.890 


P 


11 211-2 




-21.3 


-22.3 


60314 


44.59 


9'! c24.0 


1.080 


P 


10.3l0i 


49 8^^-^ 
^^•^49.5 


-19.9 


-20.8 


60421 


44.24 


22 623.1 

^^•"20.4 


1.600 


P 




^18.5 


-19.9 


-20.7 


60422 


45.17 


23 023-3 

^"^•^22.6 


1.890 


P 




71 3'7i-^ 


-21.1 


-21.6 


60423 


44.34 


99 922.8 
^^•^20.4 


1.560 


P 


io.9i«:^ 


007 q393.6 


-22.6 


-23.3 
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Table 2. Morphological parameters of Type-2 QSO hosts at z< 1.2. The morphological classes in column 4 are : 1 = bulge domi- 
nated; 2 = disk dominated; 3 = possible merger. 



XID 


z 


magi 
AB 


Class 


16 

ri/^ 
pixels 


Rp 
pixels 


e 


C 


A 


G 


M20 


41 


0.962 


21.73 


1 


3.01 


6.01 


0.10 


2.99 


0.20 


0.53 


-1.70 


70 


0.688 


20.51 


2 


15.14 


30.35 


0.28 


2.42 


0.15 


0.47 


-1.66 


151 


0.797 


21.99 


3 


21.17 


51.81 


0.42 


2.36 


0.45 


0.58 


-0.66 


200 


1.156 


23.34 


3 


11.18 


31.52 


0.25 


3.24 


0.31 


0.49 


-1.89 


211 


1.188 


23.94 


1 


4.91 


13.06 


0.01 


2.99 


0.10 


0.60 


-1.83 


212 


0.931 


22.47 


1 


7.73 


19.87 


0.17 


3.16 


0.07 


0.60 


-2.10 


229 


0.864 


21.67 


1 


8.23 


18.33 


0.20 


2.72 


0.06 


0.55 


-1.92 


292 


0.618 


21.68 


2 


21.47 


41.24 


0.67 


2.45 


0.18 


0.53 


-1.79 


313 


0.970 


22.51 


2 


10.03 


22.59 


0.40 


2.77 


0.12 


0.53 


-1.84 


323 


0.839 


21.29 


2 


16.77 


36.11 


0.56 


3.49 


0.15 


0.51 


-2.29 


411 


0.952 


22.56 


1 


6.52 


15.84 


0.23 


2.95 


0.08 


0.54 


-1.98 


413 


1.023 


21.46 


3 


31.00 


83.23 


0.31 


3.51 


0.39 


0.50 


-1.86 


2208 


1.130 


23.04 


1 


7.99 


21.86 


0.25 


3.35 


0.23 


0.58 


-1.90 


2210 


0.968 


21.48 


1 


13.38 


35.99 


0.18 


3.75 


0.10 


0.57 


-2.28 


2237 


0.944 


20.94 


2 


12.57 


34.64 


0.44 


3.97 


0.14 


0.58 


-2.33 


2289 


0.833 


20.79 


1 


6.68 


15.90 


0.29 


3.69 


0.09 


0.61 


-2.18 


2440 


1.175 


22.27 


1 


9.62 


22.94 


0.19 


3.05 


0.12 


0.58 


-1.92 


2597 


1.119 


22.89 


1 


9.60 


29.59 


0.20 


3.52 


0.20 


0.60 


-2.06 


5033 


0.670 


21.02 


3 


18.95 


39.91 


0.20 


1.99 


0.21 


0.53 


-1.17 


5053 


1.200 


23.83 


1 


8.78 


29.95 


0.29 


3.25 


0.24 


0.53 


-2.06 


5229 


1.030 


23.32 


1 


6.46 


14.65 


0.23 


2.95 


0.13 


0.56 


-1.90 


5288 


1.000 


21.98 


1 


7.52 


18.34 


0.26 


3.09 


0.06 


0.57 


-1.96 


5326 


0.950 


23.97 


3 


7.76 


17.84 


0.52 


2.69 


0.32 


0.56 


-0.99 


5400 


0.680 


21.18 


1 


10.04 


23.91 


0.20 


3.42 


0.07 


0.62 


-2.18 


5414 


0.910 


22.59 


1 


8.60 


19.75 


0.46 


2.75 


0.09 


0.54 


-1.81 


5427 


1.177 


22.13 


2 


18.08 


45.91 


0.26 


3.24 


0.21 


0.54 


-2.11 


5441 


1.170 


22.77 


1 


14.84 


38.39 


0.30 


3.81 


0.24 


0.51 


-2.16 


5467 


0.870 


23.13 


1 


6.63 


15.94 


0.11 


2.94 


0.11 


0.57 


-1.94 


5496 


0.694 


21.06 


3 


12.10 


27.63 


0.41 


3.28 


0.23 


0.59 


-1.84 


5546 


1.200 


22.25 


1 


4.16 


8.85 


0.10 


2.80 


0.14 


0.56 


-1.81 


5547 


0.960 


23.77 


2 


6.37 


20.50 


0.15 


3.75 


0.20 


0.59 


-2.14 


5548 


1.060 


22.77 


2 


13.83 


30.79 


0.48 


3.21 


0.16 


0.53 


-2.14 


54514 


0.707 


20.41 


1 


13.76 


42.68 


0.32 


3.72 


0.11 


0.66 


-2.12 


60314 


1.080 


23.26 


1 


6.51 


15.56 


0.24 


2.69 


0.13 


0.55 


-1.81 
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Appendix A: SED fits 

SED decomposition for the sample of Type-2 QSOs. Black circles are the observed photometry in the rest-frame. Purple and blue 
lines correspond respectively to the galaxy and AGN template found as the best fit solution, while the black line shows their sum. 
The triangle is the predicted AGN emission at 12.3 jim based on the Gandhi et al. (2009) correlation. We tried to anchor the AGN 
template to this value allowing an uncertainty egual to the dispersion observed by lGandhi et al. I (l2009l) in their correlation. 



